;performs stddev of 2d array to return
; along the axis

function stddev_2, imgh
n_imgh=(size(imgh, /dim))[1]
means=total(imgh, 2)/n_imgh
sqrs=total(imgh^2, 2)/n_imgh
return, sqrt(sqrs-means^2)
end
